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Abstract. We study the Casimir energy due to bulk loops ol matter fields in codimension- 
two brane worlds and discuss how effective field theory methods allow us to use this result 
to renormalize the bulk and brane operators. In the calculation we explicitly sum over the 
Kaluza-Klein (KK) states with a new convenient method, which is based on a combined use of 
zeta function and dimensional regularization. Among the general class of models we consider 
^ \ we include a supersymmetric example, 6D gauged chiral supergravity. Although much of our 

discussion is more general, we treat in some detail a class of compactifications, where the extra 
dimensions parametrize a rugby ball shaped space with size stabilized by a bulk magnetic flux. 
The rugby ball geometry requires two branes, which can host the Standard Model fields and 
carry both tension and magnetic fiux (of the bulk gauge field), the leading terms in a derivative 
expansion. The brane properties have an impact on the KK spectrum and therefore on the 
QQ , Casimir energy as well as on the renormalization of the brane operators. A very interesting 

1^-^ ■ feature is that when the two branes carry exactly the same amount of flux, one half of the 

bulk supersymmetries survives after the compactification, even if the brane tensions are large. 
We also discuss the implications of these calculations for the natural value of the cosmological 
constant when the bulk has two large extra dimensions and the bulk supersymmetry is partially 
preserved (or completely broken). 



X: 

H . 

Cd . 1. Introduction 



Higher dimensional theories have a variety of motivations in the physics of fundamental 
interactions. The original one, in the pioneer works of Kaluza and Klein (KK), was the 
unification of the known forces. This idea survived as the decades passed by through superstring 
theory, which provides a framework where gravity can be unified with the other interactions and 
strongly indicates that the space-time contains a number of extra dimensions. Other advantages 
of higher dimensional setups include the possibility to understand the structure of fermion masses 
and couplings [1] in the Standard Model. 

Remarkably extra dimensions also provide mechanisms to understand technical naturalness 
problems. For example, the gauge hierarchy problem, which queries why the Fermi scale is 
so much smaller than the Planck mass, can be addressed in the large extra dimension (LED) 
scenario or by so called warped geometries. It is then natural to ask whether extra dimensions 
can help us understand the most serious fine tuning issue we know, the cosmological constant 
problem [2]. 



To understand the origin of such problem it is useful to look at the usual Einstein-Hilbert 
gravit}(3 coupled to matter 

^ \ , ^matter \ 



V=5 2k| 

where Cmatter is the part which depends on the matter fields. This represents the leading 
gravitational tensor theory in the large wavelength limit. Here K4 is the 4D Planck scale and Aq 
the tree level value of the cosmological constant. From this equation it is clear that the quantum 
vacuum energy (p) produced by the gravitational field g^y and any matter particle contributes 
to the cosmological constant A = Aq + (p). The cosmological constant problem consists in the 
mismatch between the observed value A ~ (10~^eV)^ and the (relatively) enormous contribution 
to {p) due to the known forms of matter: any Standard Model particle with mass m adds to A 
a quantity of order rn^ and therefore a huge fine tuning of Aq is required to obtain the observed 
value. 

Any solution of such problem must involve a modification of gravity at an energy scale 
roughly of order 10~^eV. Indeed we know all non-gravitational forces up to the ~ 10 TeV scale 
and the only possible modification should therefore occur in the gravity sector, which is much 
less constrained by the experiments; if such modification emerges only at some intermediate 
energy scale Eint, between 10~^eV and 10 TeV, then the low energy effective theory at energies 
much smaller than Eint would be again the Standard Model plus Einstein's gravitational theory 
and we would be back to the original problem. 

The required modification occurs in any model with two large extra dimensions. Let 
us see why. The LED scenario predicts the following ratio between the 4D Planck mass, 
Mpi = 1/(\/2k4), and its higher dimensional counterpart, M = 1/(\/2k)^/(^~^), [3] 

# = (^o^-^ (2) 

where D is the full space-time dimension and r is the linear size of the D — 4 dimensional space 
volume (that is the length scale of the extra dimensions). In this framework one addresses the 
gauge hierarchy problem by choosing roughly M ~ TeV, which for D = 6 gives us a KK scale 
of order 1/r ~ 10~^eV (above which gravity is modified). It is important to notice that a true 
solution of the gauge hierarchy problem requires a mechanism to dynamically stabilize r to this 
large value; one way to do so is to consider a compactification on a topologically non-trivial space 
with a magnetic flux (flux stabilization). Also, notice that in order for the extra dimensions 
to be so large, only gravity (among the known interactions) can propagate in the bulk while 
the Standard Model fields should be confined on a 3-1-1 dimensional brane (3-brane henceforth), 
which in the case of interest here, D = 6, has to be of codimension-two. 

Of course, the fact that gravity is modified at the required energy scale is not by itself a 
solution of the cosmological constant problem. We need a further mechanism protecting A from 
large quantum corrections. Supersymmetry is a possible candidate because the fermion and 
boson contributions to the vacuum energy cancel exactly if supersymmetry is unbroken, or, if 
supersymmetry is broken at an energy scale mg, produce a net vacuum energy of order m|. A 
crucial observation now is that the supersymmetry breaking scale in the bulk does not have to 
be of the same order of that on the brane, which is required to be larger than the TeV scale 
by the LHC experiments. A supersymmetry breaking brane Lagrangian 5Lij, which, including 
quantum corrections, is expected to be of order TeV gives rise to a supersymmetry breaking 
splitting 5m\j^ in the KK spectrum of order K^bL\,lr^. Here the factor reminds us that 
the brane physics is communicated to the bulk through gravitational interactions while 1/r^ is 



^ We adopt the mostly plus signature and the curvature conventions of Ref . [3] . 



there for dimensional reasons. To the extent that m| is given by Smj^j^ we therefore obtain a 
cosmological constant of the correct order of magnitude. To summarize 



m,^6m^^ -2 -2 



where the question mark reminds us our expectation. To confirm it exphcit calculations in a 
concrete model are therefore needed. 

The idea that supersymmetric large extra dimensions with codimension-two branes can 
address the gauge hierarchy and the cosmological constant problem was originally proposed 
in [5]. We refer to these works for an extended discussion. But the fact that extra dimensions 
in general can help with the cosmological constant problem is older (see for example Refs. [6]) 

Since any issue of technical naturalness is a quantum problem we believe it is important to 
have a sistematic approach to compute quantum corrections in codimension-two brane worlds. 
Here we discuss quantum corrections due to bulk (massive) states, focusing on the one-loop 
approximation. The main subject of this article is the Casimir energy produced by integrating 
out a (massive) bulk field, and how to obtain it from its KK mass spectrum. Much of our 
discussions is based on Refs. [71 El EJ [TUl E], but here we also obtain some original results. 
The main one is a convenient way to compute Casimir energies from KK spectra, which is 
based on a combined use of dimensional regularization (for the ultraviolet divergences) and zeta 
function regularization (for the KK sums). The Casimir energy can then be used to calculate 
the renormalization group equations (RGEs) of both bulk and brane coefficients in the quantum 
action. 

Although some results can be applied to any codimension-two model, we study in some detail 
a concrete class of compactifications, where the extra dimensions form a rugby ball shaped space 
and the flux of a gauge field stabilizes their size, r. These solutions are supported by two 3- 
branes having both tension [12] and magnetic flux (the leading terms in a derivative expansion 
of the brane action) [13j. We discuss a class of models having these solutions, including one with 
bulk super symmetry, 6D gauged chiral supergravity [141 [T5l I16j . which we refer to throughout 
this paper as our supersymmetric example. As shown in Ref. when the two 3-branes 

have identical localized fluxes one half of the bulk supersymetries are preserved (if the bulk is 
supersymmetric), implying, interestingly, that the 4D vacuum energy vanishes in this limit. 

For rugby ball compactifications, the KK towers of many bulk fields are known from previous 
calculations [3 El [9l \TU[ |11] and we make use of these results to compute explicitly the Casimir 
energy and the renormalization of bulk and brane coefficients as a function of the bulk mass m 
of the field we integrate out. Our computational method is very convenient and confirms the 
results of Ref. [10]. The final expression for the Casimir energy is a polynomial function in mr 
of degree 6, where the coefficients depend on the brane tension and the bulk and brane localized 
fluxes. In the supersymmetric case, all these coefficients vanish for identical localized fluxes or, 
more generically, they are suppressed by the difference between these fluxes. As shown in [11] 
the 4D cosmological constant inherits this suppression and can be of the observed size. 

Let us give the outline of this article. In section[2]we introduce our family of models, including 
both the bulk and brane actions, and the rugby ball solutions; as a concrete supersymmetric 
example we define 6D gauged chiral supergravity and reexamine the technical naturalness of the 
cosmological constant in this specific case. We then review the KK spectra of scalars, fermions 
and vectors for rugby ball compactifications in section El The Casimir energy calculation and 
the bulk renormalization is performed in section HI treating in some detail the codimension-two 
case and, more specifically, the rugby ball solutions. In section Owe finally discuss the form of 
the Casimir energy due to loops of specific matter fields: the simplest case of a real scalar and 
(massive) matter multiplets of 6D supergravity. 



2. The class of models 

We focus on a class of models which include, in addition to the metric tensor Qmm, a set of gauge 
fields A'^j, scalars (p^ and fermions ip'^. The bosonic part of the Lagrangian for these fields is 

-j= = - - \ O^MDM<P' D''4^ - ^-HabicP) Ft,,F''^' - V{<P) , (3) 



where is the gauge-covariant derivatives for the scalars, F'^j^ is the field strength of 
and 'Hab{4>) aiid V{(t)) are generic functions of the scalars. 

This is general enough to describe the linearized dynamics of matter supermultiplets in 
explicit 6D supergravities. The supersymmetric example we shall refer to throughout this paper 
is 6D gauged chiral supergravity. However, the general class of models we consider also includes 
non-supersymmetric theories and the results we shall discuss in the following sections hold 
for them as well (unless otherwise stated). The supersymmetric field content is given by a 
supergravity-tensor multiplet {qmn, Bujpf, cj), -0^/, x) - with metric tensor, anti-symmetric Kalb- 
Ramond field (with field strength Gmnp), dilaton, gravitino and dilatino - coupled to some gauge 
multiplets {A'^j, A) - with gauge fields and gauginos - and some hypermultiplets ($^, - each 
of them having hyperscalars and their chiral hyperini. The fermions are all Weyl spinors and 
satisfy T^tpj^j = tpM, ^7X = ~Xi TyA = A and Fy^' = — We will consider a matter content 
with gauge group of the form Q = Q x [/(l)/j, where C/(1)_r is an Abelian R-symmetry, and G is 
a generic product of simple groups. The bosonic Lagrangian of 6D gauged chiral supergravity 
is given by 



I At* 



where g is the \J(X)r gauge coupling and we have set Gmnp = for simplicity. 

As we have already mentioned a realistic realization of the large extra dimensions idea requires 
the presence of 3-branes, one of which supporting the Standard Model fields. We take the brane 
action to be [T3] . 

Sb = - J (fx^/^Li, 
with Lb = n-^e"'^Fmn + — R + -{d^f + ---, (5) 

where gb is the gauge coupling of the gauge field corresponding to the second term in ([5]) (which 
has to be Abelian in order for that term to be gauge invariant). The ellipses denote other terms 
involving two or more derivatives (including in principle the Standard Model fields), and Tb, Ab, 
Bb, Cb and so on could depend on the scalars 

One can now ask again whether the observed value of the cosmological constant emerges 
naturally in this framework. While no positive answer is found in non-supersymmetric cases, 
the supersymmetric model has been proven to have very interesting features OUZj. One can 
compute the cosmological constant in three steps: first, the brane fields are integrated out at 
the quantum level; second, the classical integration in the bulk is performed; third, the quantum 
corrections of the bulk integration are computed. The first and second steps always give a 
vanishing cosmological constant if the dilaton does not couple directly to the branes [SJ [18] . Ref . 
|17j has recently shown that the third step is not dangerous either, at least for representations 
of the supersymmetry algebra which are massless in the 6D sense. The essence of the argument 
is that e^*^ acts as a loop counting quantity (as it can be proved by going to the frame defined 
by gMN — ^'^gMn) and is very small, of order l/(r^TeV^). The last property can be understood 



by noticing that the bulk Lagrangian enjoys a classical scale invariance under which c'^Qmn does 
not change; in the large extra dimension setup the value e'^r'^ is therefore expected to be fixed 
by the field equations to a value of order 1/TeV^. The contents of the following sections are 
useful tools to address the same question in the presence of massive 6D supermultiplets (see also 
[T9] for a study of bulk UV sensitivity for Ricci-flat geometries, including the gravity sector, but 
without branes.). 

Let us conclude this section by discussing the solutions of the models we have considered. 
The simplest solution preserving 4D Poincare invariance and providing a flux stabilization is a 
{Minkowski)^ x S"^ compactification, where the metric of the two extra dimensions is that of a 
sphere with radius r, 

ds^ = r\de^ + sin^ e dip^) . (6) 

[0 < 9 < TT and < f < 2tt). Also, the scalars 0* are constant and the non- vanishing components 
of the field strength are Fmn = f^mm where / is a constant. If there are some matter fields 
having a non-trivial chargcl, qg, under this background gauge field, then its field strength F^n 
has to satisfy the following quantization condition [20] 

IttN = q F = Anr'^qf , (without brane sources) (7) 

where = 0, ±1, ... is an arbitrary integer. Then / must satisfy 

N 

f = ^ . (without brane sources) (8) 

In the super symmetric model, the field equations also imply 

'"^ and / = ±- ^ . (supersymmetric case) (9) 



4^2 J.2 2gr^ 

Since all the parameters of the bulk theory are expected to be of order TeV (to the appropriate 
power) in the large extra dimensions setup, we notice that the first of these two conditions 
confirms our general expectation that e'^r^ is of order 1/TeV^. The second one, combined with 
Eq. ([8]), gives us = ±qg/g. When the background gauge field is along the C/(1)r, and so 
g = g, N = ±1 and this solution preserves half of the supersymmetries of the bulk theory 
|16j . For other embeddings of the background gauge field supersymmetry is instead completely 
broken. 

A simple way to introduce 3-branes in this case is starting from the sphere metric and 
demanding ip to have a modified period, (/? ~ + 27ra, where a is a positive constant [12^J. This 
procedure introduces two conical singularities, one at the north pole and the other one at the 
south pole, with identical defect angles, 5 = 27r(l — a). The resulting internal space is called 
the rugby ball. 

The brane action in ([5]) can support these conical singularities and the near-brane boundary 
conditions imply [13J that the defect angle is 

Sb = K^Lb . (10) 

In order for this angle to be positive we will assume that Lf, > at the background solution. 
The brane sources also changes the flux quantization condition. This arises because the branes 
themselves can support a localized flux: 



^ We denote with g the gauge couphng of the background gauge field, while q is the charge of a given field, which 
we assume to be integer. 



If the branes are identical the total flux localized in this way is ^ := J2b ^b, in terms of which 
the flux-quantization condition is (when the background gauge fleld has some localized flux) 



2tt^ + ATTar'^qf, 



(12) 



where 5*2 (a) is the rugby ball. 



The normalization constant / is then given by 



/ 



(with brane sources) 



(13) 
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where M := uj{N — $) and u := 1/a. 

In the super symmetric case, the above-mentioned classical scale invariance can be broken by 
the boundary localized fluxes, the second term in Eq. ([5]), which gives mass to the corresponding 
flat direction |13] : another reason why these fluxes are useful is that they allow us to recover 
supersymmetry, when the background gauge fleld is along the U{1)r, in a continuous limit 
[131 [T7]. To understand this let us take the rugby ball solution, g = g and <^fe = 0, then = ojN 
and the second condition in Eq. ([9]) together with ()13p implies uoN = itg, which does not allow 
you to approach the super symmetric value w = 1 continuously. In the absence of instead we 
do not find this obstruction: using the second condition of ([9]) in ()13p this time we obtain 



In the presence of boundary localized fluxes supersymmetry can be broken by an arbitrarily 
small amount. For this reason in the super symmetric model we will always take the boundary 
localized flux in the f7(l)_R direction, gi, = g. 

Let us conclude this section by mentioning that, for branes carrying both tension and flux, 
there is, remarkably, a case in which one half of the 6D supersymmetries is preserved; this occurs 
when the background gauge fleld is along the U{1)r generator, for which there is exactly the 
same amount of localized flux on the two branes [llj (i.e. = $/2 both at the north and south 
branes). Such property has important implications regarding the smallness of the cosmological 
constant. 

3. Spectrum in codimension-two brane worlds 

We now move to the analysis of the linear perturbations around the 3-brane solutions we 
have considered. Since these configurations preserve 4D Poncare symmetry, such analysis is 
equivalent to computing the 4D particle spectrum (defined in the usual sense). This will provide 
us with additional physical information; it is an important computation to study the stability 
of the background solutions |8t i9j and, of special relevance for this article, it is an intermediate 
step to determine the quantum corrections to the 4D vacuum energy. Indeed in dimensional 
regularization the one-loop contribution to the quantum potential (the Casimir energy) due to 
a generic fleld with bulk mass m is 




(14) 



which, in the case g = g and N/q = ±1, gives 




(15) 




(16) 



where the coUective index n includes all KK numbers and m„ represents the full set of KK 
masses. Also (-1)^ is 1 for bosons and —1 for fermions. The 4D particle spectrum is therefore 
an important ingredient to compute V and we will use Eq. (|16p to compute the renormalization 
of bulk and brane couplings in sections H] and [5j 

In the presence of two extra dimensions we have two KK numbers, n = j, n; in the rest of 
this section we give the form of mj^ := Xjn/i''^ for (minimally coupled) scalars, fermions and 
gauge fields on top of the rugby ball geometry sourced by branes with both tension T^, and flux 
$f, (the leading terms in a derivative expansion). In the absence of localized fluxes, these 
were compute 41 in Refs. [71 El IB]. Subsequently, theh for m in the presence of was derived 
in do]. 

Let us consider flrst the simple case of a minimally coupled real scalar, satisfying the equation 
{DniD^' + m'^)(j) = 0, that is coupled to the background gauge field through a monopole number 

and boundary localized fiuxes ^f,. In this case the scalar spectrum (in the north patch of the 
gauge potential) is 



jn 



7 H n 
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2 ' 



(17) 



where 



i = 0,l,2,..., n = 0,±l,±2 



and {^-) is the fiux localized on the north (south) pole of the rugby ball, where cos0 = +1 
(-!)• 

Moving to fermions, the KK spectrum for a field, satisfying {Ij) + m)il) = 0, that is charged 
under the U{1) is (using again the north patch of the gauge potential) 



A 



'jn 



UJ 

+ 2 



"-1/2 



7V + $_ + — 
2uj 



1\2 



where ni/2 = n — a/2 and a G {±1} corresponds to the 4D helicity of the spinor, of which there 
are 2 (4) each for a 6D Weyl (Dirac) spinor. 

Finally let us give the KK spectrum arising from a gauge field, which assume not to be in 
the Lie algebra direction where the background gauge field hefl There a two cases: the gauge 
field can either be massless or massive. We begin with the massless case, when the field satisfy 
the equation g'^^ DmF^p = 0. In an appropriate gauge (e.g. light-cone gauge [25t [26 l [HI |9]) 
the 6D gauge field can be decomposed into four components, each with a spectrum (once again 
in the north patch of the gauge potential) given by 



n 



+ 2 



n - iV + - - 



1 

+ 2 



(l+A/'2) 



4 



(19) 



where G {0,0, +1,-1} for each of the four components and we assume N = 0, ±1 to ensure 
stability [271 E]. We observe that two modes have exactly the same spectrum as scalars (i.e. those 
with 1^ = 0), while the other two (with ^ = ±1) have almost the same spectrum. For massive 
gauge fields we have a more complicated situation because we need a scalar field ^ that is charged 



^ Those works considered a more general compactification [JJ than the rugby ball one, which involves generically 
non-factorizable geometries. Here, however, we focus on the rugby ball case for the sake of simplicity (see [22] for 
the explicit unwarped limit). 

* This assumption simplifies our calculation because the gauge field perturbations along the same Lie algebra 
direction as the background gauge field mix with the metric perturbations. See, however, Ref. _9!, for a discussion 
of this more complicated sector, but in the absence of boundary localized fluxes. 



under the gauge field we are studying in order to give mass through the Higgs mechanism. 
In order to interpret (<&) 7^ as a 6D spontaneous symmetry breaking, we require (<I>) to be 
constant and to be at the minimum of U. Then, in order to solve the background scalar equation, 
D^D^'^ = 0, we also demand that ($) 7^ does not break the U{1) where the background 
gauge field lies: otherwise it would not be possible to have (<&) constant, at least in the sphere 
compactification of interest in this paper [281 EHl ES]. If we choose again the light-cone gauge, 
it is possible to show that a massive gauge field leads to the 4D spectrum of a massless gauge 
field, Eq. ([H]), plus that of a scalar, Eq. (fT7|) . 



4. Casimir energy and bulk renormalization 

The purpose of this section is to compute the Casimir energy due to bulk loops and to show how, 
consequently, the bulk and brane couplings renormalize. This is an important step to address 
any problem of technical naturalness, such as the gauge hierarchy and the cosmological constant 
problem. We will consider in some detail the case in which the branes are of codimension-two, 
but some of our results will be valid in more general brane worlds. The renormalized couplings 
will depend as usual on a renormalization energy // and we will compute explicitly their RGEs. 



^.1. A general technique to compute the Casimir energy 

The starting point of this calculation is the formula for the one-loop contribution of a single 
field to the quantum potential in Eq. ()16p . Notice that, modulo terms which are independent 
of X, we have InX = — J^(ds/s) exp(— sX) and therefore 



y=-l /-'(-I)" i:J^.ri -P [- (/ + + -^)] . (20) 

Performing the integral in d'^p and rescaling the variable s we obtain 
V 



Consider first the case Re(d) < and rewrite the integral in Eq. (|21|) as 

We now extend this integral function by analytic continuation to all complex d except the 
non-positive integers (where the function has simple poles). By using 

r{-d/2) = -j^ + l-^ + Oid-4) , (23) 
where 7 is Euler's constant, one then obtains 

V = - ^ ^2^, {-If E + ^'^y^' (- + + -•) ' (24) 

where the dots represent finite r-independent terms. 

One possible renormalization scheme (which we will adopt from now on) is to subtract the 
divergent part in the brackets of Eq. (j24p . The renormalized potential Vr (the Casimir energy) 
can then be written as 



where Tq has to be identified as a ultraviolet (UV) scale which can be computed once the UV 
completion is known and 

C:=-\{-lYY.i^'ml + rVf\ (26) 

n 

It is important to notice that all we need in order to compute this coefficient is the divergent 
part of y, as it is clear from Eq. ([2l|) . 

The quantity on the right hand side of Eq. ()26p is divergent and has to be regularized. 
Notice that the exponent d/2 can effectively act as a regulator for the sum over n (zeta function 
regularization [30J): for d = 4 the sum is divergent, but one can (and we will) compute (|26|) for 
those d such that the sum is convergent and then consider the analytic continuation of the final 
result at d = 4. 



^.2. The Casimir energy for rugby ball compactifications 

A generic form of the spectrum which covers all the cases encountered here, Eqs. ()17p . (jlSp and 
(USD, is 



m 



'jn ■— ^2 '^i" ' 



where A 
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J + — |n + 6+1 + — |n - 6_| + a 



(27) 



and b±, a and r are real parameters; the only assumption we make is that b± and a are 
independent of j and r is independent of both j and n. 

Then the contribution to the C parameter from a single field is 



C 



-11 



^ "I [i + ^ |n + 6+1 + ^ |n - 6_| + a J - r + (mr) 



d/2 



(28) 



where (—1)^ = 1 for a boson and (—1)^ = —1 for a fermion. We can now expand {...}^^^ by 
using the binomial series to obtain 



C 



-ly ^ T{k-d/2) 

2 t^.t^o^'r(-d/2) 



J + — |n + 6+1 + — |n - 6_| + a 



d-2k 



T — [mrY 



(29) 



The sum over j can be performed by means of the following representation of the Hurwitz zeta 
function C(S)C) (valid for Re(s) < and Re(c) > 0) 



C(-s,c) :=^(j + c) 



1 



i=o 



T[-s] 







oo y-l-s Q—cy 
dy . 

^ 1-e-y 



(30) 



Indeed setting s = d — 2k and c = c„ = (|n + 6+| + \n — 6_|)a;/2 + a in ([30|) we have 



C 



i-lf g T{k - d/2) [t - (mr) 



2^'^ roo y-l+2k-d 



2 ^ k\T{-d/2)r{2k - d) Jo ^ 1 - e^y 



dy 



(31) 



Observing that integrals of the form /q°° f_^^-y can be computed explicitly for Re(/i) > 0, 



dy 



ye 
1 - e-y 



r[l + /i]Lii+,,(l), Lii+;,(1) := E Im ' 



(32) 



fc=i 



we Taylor-expand ye^ Xln ^ "^"^ around y = 0, 

,e^^e--^=E^/> (33) 

n k'=0 • 

and obtain 

^ — kT{-d/m2k - d) -^^^-^ -d + 2k + k )Li_i_,+2.+.'(l). (34) 

We now need to take the limit d — t- 4. The previous expression turns out to be well defined in 
this limit. Also, because of the T(—d/2) in the denominator, only a finite number of k and k' 
contributes to the sum: A; = 0, 1, 2, 3 and k' = 0, 1, 2, 3, 4, 5, 6. This implies two things: (i) C is 
computable once we know X^n^"'^"^' (^0 ^ generic structure 

C = ^-^{mr)^ — -^(mr)^ + si(mr)^ — S2 ■ (35) 

Since the coefficient C can be computed from the mere knowledge of the UV divergent 
part of the Casimir energy, we should expect that this formula can be applied to codimension- 
two compactifications which are more general than the rugby ball one. This is because UV 
divergences are related to the local structure of the space-time and are therefore insensitive to 
global properties such as its topology. 



4-3. Bulk and brane counterterms and renormalized couplings 

The models we have considered are non-renormalizable and so not all divergences can be 
reabsorbed in counterterms of the same form as the terms in the classical action. However, 
one can show in very general terms that the number of counterterms needed is always finite at a 
given order in perturbation theory. In this subsection we review [IQ] the renormalization due to 
bulk loops in codimension-two braneworlds (and in particular for rugby ball compactifications) 
at the one-loop level. At the end we will therefore obtain a finite number of counterterms and 
renormalized couplings. To perform this calculation we will use effective field theory methods 
(see for example [31]). 

Although only bulk loop^ are computed here, both bulk ad brane counterterms are needed. 
There is an important difference between them. The bulk counterterms, unlike the brane ones, 
do not depend on the brane properties and so they can be computed in the sphere limit, a — )• 1. 



4-3.1. Renormalization of the bulk interactions. To capture all the terms needed to reabsorb the 
UV divergences we write down the most general local Lagrangian with the chosen field content 
and set of symmetries, which we organize in a derivative expansion C^g = £^*q -|- £^*2 + ^b4- 
After renormalization this generates a corresponding series of renormalized interactions £^ = 
C^gQ + CJ'g2 + ^i34- Focusing on the fields that are non-zero in the background we have 
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^ Loops of brane localized fields can be calculated and the corresponding divergences can be reabsorbed with the 
usual 4D field theory methods. 



where B? {B?) is a generic linear combination of terms which are quadratic (cubic) in the 
curvature, that is 



o-R R + 2b R RmnR'^'^ + Cn RmnpqR'^"^^'^ (40) 

where Qr + + Cr = 1 so that R^ = R? when speciahzed to the sphere geometry (for which 
RmnpqR^^^'^ = 2RmnR^^ = R^ = 4/r*^). A similar expression is used for . Calculations on a 
sphere can only provide the overall couplings , and not the separate parameters such as 
ajj, bji and Cr (see however [10] and references therein to know the latter quantities). 

Evaluating the renormalized action at the background sphere solution and integrating over 
the extra dimensions gives 
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Therefore A, Cr, Cr'^j and C^a can be read off respectively from the r^, r^, r~^, and terms 
in Vr (see Eqs. ([25]) and ([55]) ). while the Ca and Car coefficients are identified as the J\f'^/r'^ 
and Af^/r^ terms respectively. This implies that integrating a bulk field with mass m gives the 
following contribution to the RGEs 
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(44) 

{A^fr^P '2 - (47r)W2 

and so on. The "sph" in s^^^' ^ emphasizes that these quantities are evaluated on the sphere, while 
the superscript '/c' denotes terms involving k powers of N . The renormalization of the gauge- 
field terms, Ca and Car-, has been done by looking at the AA-dependent divergences produced 
when a particle with charge qg runs in the loop. 
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4-3.2. Renormalization of the brane interactions. In this case we have a dependence on the 
boundary conditions used near the brane but the result should be independent of the boundary 
conditions on distant branes. 

To compute the brane contributions we first subtract the (boundary condition independent) 
bulk contributions found above. Noticing that the bulk counterterms should be integrated over 
the volume of the rugby ball, which is Airar'^, we define 

<5sr = s,-as^P\ (46) 

and use 6s\"^ = to extract how the interactions on each individual brane renormalize. 

This can be done as before, by distinguishing the interactions that depend on the gauge field 
which is non-zero on the background from those that do not. In the following we understand 
the label (6) in to have a simpler notation. 



Writing the most general local brane Lagrangian organized in a derivative expansion, 
C\ = C^Q + + + + •••! and dropping terms that vanish at the background, we have 



C 



bi 



-62 



-7 



Cxb mn p 



K 4(7^ 



63 



-7 



64 



-7 



2^2 

Cjj26 + 



.4ii 6 r> ,mn 77 



2 , '^^ Caa 6 



r? 17 TpMN 



(47) 
(48) 
(49) 
(50) 

(51) 



and so on, where '^^y := g^i^d^x'^' d^x^ (with the right-hand side computed at the brane 
position) is the induced metric on the brane. 

Evaluating these at the background solution gives the following contribution to the Casimir 
energy 
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Using Eqs. (|25p and ()35p then gives the following RGEs 
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where Js^ are terms with k powers of A/". 
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5. Casimir energy for rugby balls in explicit cases 

In this section we apply the method of sections 14.11 and 14.21 to compute the Casimir energy, i.e. 
the coefficients Sj, produced by specific matter fields, for rugby ball compactifications sourced 
by branes with tension and fiux. Indeed the Sj is all we need to obtain the renormalization of 
the bulk ad brane couplings, as it is clear from Eqs. ()42p - (|45p and ()53p . Moreover, as we shall 
comment later on, the Si can be used to extract the 4D cosmological constant [10\ lllj. 



5.1. A single real scalar. 

We observe that the method of section S2] can be applied in the case of real scalars because the 
spectrum in (|17p has the form ()27p with 

1 1+7^2 

6+ = -<&+, 5_ = iV-<&_, a = -, T = ^^. (54) 



The (renormalized) Casimir energy produced by a real scalar is given by Eq. ()25p with C given 
in Eq. ()35p and one obtains the following Sj coefficients: 
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The method of sections 14.11 and 14.21 can also be apphed to fermions and gauge fields as their 
spectra, Eqs. ()18p and ()19p . are both of the form given in (|27p . The explicit expressions for the 
Si coefficients for fermions and (massive) gauge fields can be found in |10] . 

5.2. Supermultiplets 

Let us now consider supermultiplets of 6D gauged chiral supergravity focusing on the 
hypermultiplets and the gauge multiplets (for which we further restrict to the case in which 
the gauge field is zero in the background). The main reason is that the cancellation of gauge 
and gravitational anomalies typically require hundreds of such supermultiplets [S51 [33] and so 
their contribution is expected to dominate the Casimir energy. For example, the first found 
anomaly free theory of this sort has a large {Eq x Ej xU{l)ji) gauge symmetry with many (456) 
hypermultiplets [32} [33], 

A massless hypermultiplet consists of four massless scalars (called hyperscalars) and one 6D 
Weyl fermion, the hyperino. A massless gauge multiplets is made of one gauge field and a 6D 
Weyl fermion, the gaugino. By contrast, a massive 6D matter multiplet consists of a massive 
gauge field, a massive Dirac fermion and three scalars, a total of eight bosonic and eight fermionic 
states. Since this is also the number of degrees of freedom of a gauge plus a hypermultiplet, 
one expects to form a massive supermultiplet by having the gauge boson from a gauge multiplet 
'eat' one of the hyperscalars through the Higgs mechanism. 



5.2.1. N on- super symmetric embedding of the background gauge field. We first consider the case 
in which the background gauge field is not embedded in the U{\)ii. In this case supersymmetry 
is broken both by the branes and the bulk solution. Since the gauge field whose flux is localized 
on the branes is the C/(1)_r gauge field, the spectrum and the Casimir energy as well as the 
renormalization will not depend on <I>fe for this choice of the gauge field embedding. 



The contribution of a hypermultiplet to the Si coefficients is obtained by summing the result 
for a 6D Weyl fermion to that produced by four hyperscalar^. We obtain 
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where here N is the common monopole number of the hyperscalars and hyperino. For a massless 
[m = 0) hypermultiplet the only coefficient which matters is and the corresponding Casimir 
energy is — S2™ln(r/ro)/(47rr^)^. 

As far as the gauge multiplet is concerned, we should sum the contribution of a gauge field 
to that of a 6D Weyl fermion, obtaining for the Sj coefficients 
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where now N is the common monopole number of the gauge field and gaugino. For a massless 
gauge multiplet the only important coefficient is sf™ and the corresponding Casimir energy is 
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Therefore, if one has massless supermultiplets only, the total contribution to the Casimir 
energy is approximately given by — (s2™ + sf™) ln(r/ro)/(47rr^)^. The overall sign depends on 
the particular anomaly free model [32] that one chooses. 

For a massive multiplet made of a hypermultiplet and a gauge multiplet we have that the s. 
coefficients are sf™ 
given before, we obtain 



+ sf™. Therefore, by using the explicit expressions for s^™ and sf™ 
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From this result, and from Eqs. (j35p and ()25p . we note that for mr at most of order 1 the 
obtained Casimir energy is ~ ln(r/ro)/(47rr^)^. However, for mr ^ 1, integrating out a massive 
multiplet gives a dangerously large contribution. We will see that an extra suppression can be 
obtained when the gauge background is along the R-symmetry generator. 



The hyerscalar contribution is given by the scalar result after substituting + l/(2r'^) 'IV.. 



5. 2. 2. Super symmetric embedding of the background gauge field. Let us now turn to the case in 
which the background gauge field is along the U{1)r. Hyper and gauge multiplet contributions 
to the Casimir energy have been computed in [TT] and are quite involved. We therefore refer to 
this work for explicit expressions associated with hyper, gauge and massive multiplets. 

One important point we want to emphasize here is that s^^^ = for this background gauge 
field embedding and therefore the renormalization group equation of the bulk Casimir energy 
vanishes, iidfjVg = 0. The reason is that, as we mentioned, one half of bulk super symmetry is 
not broken in this case. 

In the particular case of identical boundary localized fluxes {^h = $/2) also the branes 
preserve one half of the 6D supersymmetries (see the discussion at the end of section [2|) . This 
implies that for balanced fluxes on the two branes also the brane Casimir energy does not run, 
^id^V^ = 0. 

When the boundary localized fluxes are unbalanced a non-trivial result arises. However, 
by continuity the result should be suppressed by the diff'erence of the two localized fluxes, 
A$. Also, as a remnant of 6D supersymmetry the coefficient s_i for massive supermultiplets 
vanishes. When the bulk mass is such that mr is at most of order one the Casimir energy is 
~ ln(r/ro)/(47rr^)^. When mr ^ 1 the contribution of a massive supermultiplet to the Casimir 
energy is of order A$(mr)'^ln(r/ro)/(47rr^)^ and can be again as small as ~ ln(r/ro)/(47rr^)^ 
for A<I> appropriately small. It is very interesting to notice that the localized fiux difference does 
not receive quantum corrections from loops involving brane localized fields only (when they are 
not charged under the corresponding gauge symmetry) and therefore taking this difference to 
be very small does not need to have the usual fine tuning. 

Let us conclude this section by mentioning that the Casimir energy is not exactly equal 
to the 4D cosmological constant, but is rather equal to the quantum action evaluated at the 
classical solution. The backreaction of the brane on the bulk is important in this case and leads 
to a sizeble correction of the solution. This point is extensively discussed in [llj . After taking 
into account this effect, however, the prediction for the (most UV sensitive part of the) 4D 
cosmological constant can remain as small as ln(r/ro)/(47rr^)^. 

6. Conclusions and outlook 

We discussed the Casimir energy as well as the renormalization of bulk and brane coefficients in 
the quantum action produced by integrating out (massive) bulk matter in codimension-two brane 
worlds. In the calculation we focused on the one-loop approximation and explicitly summed over 
the KK towers. Much of what we presented here is a review of [3 [51 El UHl E] , but we also 
provided some new results, in particular in the technique to perform the KK sums. 

Regarding the motivations, as discussed extensively in the introduction, codimension-two 
brane worlds may provide a framework to solve the gauge hierarchy and, when the bulk is 
supersymmetric, the cosmological constant problem. These are technical naturalness problems 
and as such they require sistematic ways of computing quantum corrections. 

We considered in some detail the rugby ball compactiflcations in which the size r of the extra 
dimensions is stabilized by the flux of a bulk fleld. A class of model having these conflgurations 
as solutions has been presented, including a concrete supersymmetric theory, 6D gauged chiral 
supergravity. The two 3-branes required to support these solutions can not only carry tension, 
but also a localized ffux of the same gauge field which stabilizes the extra dimensions (tensions 
and localized fiuxes are the leading terms in a derivative expansion of the brane Lagrangians, 
Eq. dSD). 

As discussed in [11], when the localized fluxes on the two 3-branes are identical one half of 
the bulk supersymmetries is unbroken, which implies that the Casimir energy and, remarkably, 
the 4D cosmological constant vanish. 

For the rugby ball compactiflcations the KK spectra of many types of bulk fields are known 



13 El El [TOl [TT] and we reviewed their structure for scalars, fermions and gauge fields (see Eqs. 
(fTT]) . (fTSjl and (fT9]l rescpectively). The exphcit form of the spectra allows us to have explicit 
expressions for the Casimir energy and, consequently, for the renormalization group equations, 
which depend on the flux and brane tensions and fluxes. The calculation of the Casimir energy 
starting from the KK spectrum exploits a novel efficient technique (see sections 14.11 and 14. 2p 
which allows us to confirm the results of |10j . 

The explicit form for the Casimir energy we obtain is a polynomial function of mr of degree 
six, see eqs. ([25]) and ([55]) . The coefficients Sj of the polynomial are computed explicitly as a 
function of the brane tensions and bulk and brane fluxes; for example, in the simple case of a 
bulk scalar they are given in Eqs. 

When the bulk particle that is integrated out is massless, or m is at most of order 1/r, the 
4D cosmological constant has the desired order of magnitude regardless of the fact that there is 
bulk supersymmetry. We show that for large m the final result can be appropriately suppressed 
if we select the super symmetric model, the bulk solution preserves one of the supersymetries 
and localized fiuxes have very similar values, such that we are close to a super symmetric setup. 

Let us mention some outlook of the results presented here. 

A possible extension of our work is the inclusion of warped geometries, which represent 
the most general solutions with 4D maximal symmetry. A first step towards this goal is the 
codimension-one case, which would be interesting by itself as the Randall-Sundrum model can 
address the gauge hierarchy problem [35] (bulk fields in the Randall-Sundrum model have been 
considered in [36]). Through the AdS/CFT correspondence quantum loops in the bulk would 
correspond to 1/Nc corrections in the CFT side, where Nc counts the number of "colors". 

Regarding again holography, we notice that the formalism of Refs. [3 [H O [ini E] to compute 
the spectrum for codimension-two brane worlds that we reviewed here can also be used to analyze 
spectral properties of holographic models in the confined phase [37]: these are obtained from 
models with one extra dimension (the holographic coordinate) by an additional compactified 
dimension and scalar, vector and fermion fields in the bulk can have a variety of uses ranging 
from condensed matteiQ [39] to quantum chromodynamics [40]. 

Finally, an interesting property of supersymmetry with two large extra dimensions is that it 
could provide a link between the observed value of the cosmological constant and the scale at 
which modifications of gravity should occur, 1/r [51]. It would be interesting to know how gravity 
gets modified in the concrete supersymmetric model we discussed. In the absence of localized 
fiuxes graviton contributions have been computed in [32], but the role of these fluxes, which are 
important for the dilaton stabilization, remains an interesting target for future research. 

A cknow ledgements 

We would like to thank Cliff Burgess, Leo van Nierop, Susha Parameswaran and Matt Williams 
for collaborations, Hyun-Min Lee for much help trying to diagonalize the supergravity sector 
in early stages of this work and Riccardo Barbieri, Oriol Pujolas, Seifallah Randjbar-Daemi 
and George Thompson for useful discussions. This work was partly supported by the EU ITN 
"Uniflcation in the LHC Era" , contract PITN-GA-2009-237920 (UNILHC) and by MIUR under 
contract 2006022501. 

References 

[1] Grossman Y and Neubert M 2000 Phys. Lett. B 474. 361 (Prepnnt |hep-ph/9912408p . Huber S J and Shaft Q 
2001 Phys. Lett. B 498 256 ( Preprint | h"ep-ph /OOlOiOSj l . Salvio A and Shaposhnikov M 2007 JHEP 0711 
037 {Preprint arXiv:0707.2455 [hep-th]). 

[2] Weinberg S 1989 Rev. Mod. Phys. 61 1. 



^ See, however, [35] for an alternative to compactification in this context. 



[3] Weinberg S 1972 Gravttation and cosmology: principles and applications of the general theory of relativity 
(New York London Sydney Toronto: John Wiley & Sons). 

[4] Arkani-Hamed N, Dimopoulos S and Dvali G R 1998 Phys. Lett. B 429 263 [Preprint h ep-ph/9803315| ). 
Antoniadis I, Arkani-Hamed N, Dimopoulos S and Dvali G R 1998 Phys. Lett. B 436 257 [Preprint 
,liep-pli/9804398). Arkani-Hamed N, Dimopoulos S and Dvali G R 1999 Phys. Rev. D 59 086004 [Preprint 
lhep-ph/9807344). 

[5] Aghababaie Y, Burgess C P, Parameswaran S L and Quevedo F 2004 Nucl. Phys. B 680 389 [Preprint 
hep-tli/0304256). Burgess C P 2004 Annals Phys. 3 13 283 [Preprint h ep-tIi/0402200 ). Burgess C P 2005 
^AIP Conf. Proc. 743 417 ( Preprint [hep-th/0 41 1 140 1) . 

[6] Rubakov V A and Shaposhnikov M E 1983 Phys. Lett. B 125 139. 

Randjbar-Daemi S and Wetterich C 1986 "Kaluza-Klein Solutions With Noncompact Internal Spaces," 
Phys. Lett. B 166 65. Gherghetta T and Shaposhnikov M E 2000 "Localizing gravity on a string - like 
defect in six-dimensions," Phys. Rev. Lett. 85, 240 [Preprint hep-th/0004014 ). 

[7] Parameswaran S L, Randjbar-Daemi S and Salvio A 2007 Nucl. Phys. B 767 54 [Preprint hep-th/0608074). 

[8] Parameswaran S L, Randjbar-Daemi S and Salvio A 2007 JHEP 0801 051 (Prepnni arXiv:0706. 1893 [hep-th]). 

[9] Parameswaran S L, Randjbar-Daemi S and Salvio A 2009 JHEP 0903 136 [Preprint arXiv:0902.0375 [hep-th]). 

[10] WiUiams M, Burgess C P, van Nierop L and Salvio A 2012 Preprint laFXiv:1210.37 53 [hep-th]. 

[11] Burgess C P, van Nierop L, Parameswaran S, Salvio A and Williams M Preprint WXiv : 1 2 10 . 5405 [hep-th] . 

[12] Carroll S M and Guica M M 2003 Preprint hep-th/0302067 Navarro I 2003 JCAP 0309 004 [Preprint 
[hep-th/0302129). Navarro I 2003 Class. Quant. Grav. 20 3603 [Preprint hep-th/0305014). Nilles H P, 
Papazoglou A and Tasinato G 2004 Nucl. Phys. B 677 405 [Preprint hep-th/0309042). Bostock P, Gregory 
R, Navarro I and Santiago J 2004 Phys. Rev. Lett. 92 221601 [Preprint hep-th/0311074). Graesser M L, 
Kile J E and Wang P 2004 Phys. Rev. D 70 024008 [Preprint hep-th/0403074). Vinet J and Cline J M 
2004 Phys. Rev. D 70 083514 [Preprint hep-th/0406141 ). Navarro I and Santiago J 2005 JHEP 0502 007 
[Preprint hep-th/0411250). Kofinas G Preprint hep-th/0506035| 

[13] Burgess C P and van Nierop L 2011 JHEP 1102 094 ( Prepnnf iarXiv: 10 12 . 26"38l [hep-th] ) . Burgess C P and 
van Nierop L 2011 JHEP 1104 078 [Preprint arXiv:1101.0152| [hep-th]). 

[14] Marcus N and Schwarz J H 1982 Phys. Lett. 115B 111. 

[15] Nishino H and Sezgin E 1984 Phys. Lett. 144B 187. 

[16] Salam A and Sezgin E 1984 Phys. Lett. 147B 47. 

[17] Burgess C P and van Nierop L 2011 [Preprint arXiv:1108.0345 [hep-th]). 

[18] Aghababaie Y, Burgess C P, Cline J M, Firouzjahi H, Parameswaran S L, Quevedo F, Tasinato G and Zavala 

I 2003 JHEP 0309 037 [Preprint ;hep-th/0308064 ). 
[19] Burgess C P and Hoover D 2007 Nucl. Phys. B 772 175 [Preprint [hep-th/0504004| ). Hoover D and Burgess C 

P 2006 JHEP 0601 058 [Preprint hep-th/0507293 ). Ghilencea D M, Hoover D, Burgess C P and Quevedo 

F 2005 JHEP 050Q 050 [Preprint hep-th/0506164 ). 

[20] Randjbar-Daemi S, Salam A and Strathdee J A 1983 Nucl. Phys. B 214 491. 

[21] Gibbons G W, Gueven R and Pope C N 2004 Phys. Lett. B 595 498 [Preprint hep-th/0307238). 
[22] Parameswaran S L, Randjbar-Daemi S and Salvio A 2010 Preprint arXiv:100 1.3271 [hep-th]. 
[23] Salvio A 2007 author's Ph.D. thesis Preprint hep-th/0701020 
[24] Randjbar-Daemi S, Salam A and Strathdee J A 1984 Nuovo Cim. B 84 167. 

[25] Randjbar-Daemi S and Sarmadi M H 1985 Phys. Lett. B 151 343. 

[26] Randjbar-Daemi S and Shaposhnikov M 2002 Nucl. Phys. B 645 188 (Prepnni hep-th/0206016 ). 

[27] Randjbar-Daemi S, Salam A and Strathdee J A 1983 Phys. Lett. B 124, 345 [Erratum-ibid. B 144, 

4 55 (1984)]. Dvah G R , Randjbar-Daemi S and Tabbash R 2002 Phys. Rev. D 65 064021 [Preprint 

arXiv:hep-ph /0102307" ) . 

[28] Randjbar-Daemi S, Salvio A and Shaposhnikov M 2006 Nucl. Phys. B 741 236 [Preprint hep-th/0601066) ; 
[29] Salvio A 2006 AIP Conf Proc. 881 58 [Preprint |hep-th/0609050p . 
[30] Hawking S W 1977 Commun. Math. Phys. 55 133^ 

[31] Burgess C P 2004 Living Rev. R el. 7 5 [Preprint [gr-qc/0311082] ). Burgess CP 2007 Ann. Rev. Nucl. Part. 

Sci. 57 329 [Preprint hep-th/0701053 ). 
[32] Randjbar-Daemi S, Salam A, Sezgin E and Strathdee J A 1985 Phys. Lett. B 151 351. 

[33] Avramis S D, Kehagias A and Randjbar-Daemi S 2005 J_ff£'P0505 057 [Preprint [hep-th/0504033 '). Avramis 
S D, Kehagias A 2005 JHEP 0510 052 [Preprint hep-th/0508172 ). 

[34] Randjbar-Daemi S and Sezgin E 2004 Nucl. Phys. B 692 346 [Preprint h ep-th/04022l7] ). 

[35] Randall L and Sundrum R 1999 Phys. Rev. Lett. 83 3370 [Preprint ;hep-ph/9905221 ). 

[36] Goldberger W D and Wise M B 1999 Phys. Rev. D 60 107505 [Preprint arXiv:hep-ph/9907218 ). 

Davoudiasl H, Hewett J L and Rizzo T G 2000 Phys. Lett. B 473 43 [Preprint arXiv:hep-ph/ 991 1262] ) . 
Gherghetta T and Pomarol A 2000 Nucl. Phys. B 586 141 [Preprint tarXiv:hep-ph/0003129| . 



[37] Witten E 1998 Adv. Theor. Math. Phys. 2 505-532. {Preprint hep-th/9803131| ). Horowitz G T, Myers R C 

1998 Phys. Rev. D59 026005. {Preprint hep-th/9808079 ). 
[38] Salvio A 2012 JHEP 1209 134 {Preprint arXiv: 1207.3800 [hep-th]]). 

[39] Nishioka T, Ryu S and Takayanagi T 2010 JHEP 1003 131 {Preprint [arXiv:091 1 .'09621 [hep-th]]). Montull M, 
Pujolas O, Salvio A and Silva P J 2011 Phys. Rev. Lett. 107 181601 {Preprint ;arXiv:1105.5392 [hep-th]]). 
Montull M, Pujolas O, Salvio A and Silva P J 2012 .JHEP 1204 135 {Preprint arXiv: 1202.0006 [hep-th]]). 
Bhattacharya J, Ogawa N, Takayanagi T and Ugajin T 2012 JHEP 1202 137 {Preprint arXiv:120 1.0?B4l 
[hep-th]]). 

[40] Basu P, Nogueira F, Rozali M, Stang J B and Van Raamsdonk M 2011 New J. Phys. 13 055001 {Preprint 

arXiv:1101.4042 [hep-th]]). 

[41] Callin P and Burgess C P 2006 Nucl. Phys. B 752 60 {Preprint [hep-ph/05l"T216] ). 
[42] Salvio A 2009 Phys. Lett. B 681 166 {Preprint [arXiv: 0909 .00231 [hep-th]]l. 



